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Transverse momentum spectra of neutral pions in the range 1 < pr < 10 GeV/c have been 
measured at mid-rapidity by the PHENIX experiment at RHIC in Au+Au collisions at ^s NN — 
200 GeV. The tv° multiplicity in central reactions is significantly below the yields measured at the 
same y/s NN in peripheral Au+Au and p + p reactions scaled by the number of nucleon-nucleon 
collisions. For the most central bin, the suppression factor is ~2.5 at pr = 2 GeV/c and increases 
to ~4-5 at pr ~ 4 GeV/c. At larger pr, the suppression remains constant within errors. The deficit 
is already apparent in semi-peripheral reactions and increases smoothly with centrality. 
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High energy collisions of heavy ions provide the means 
to study Quantum Chromodynamics (QCD) at energy 
densities where lattice calculations [1] predict a transi- 



tion from hadronic matter to a deconfined, chirally sym- 
metric plasma of quarks and gluons (QGP). The large 
center-of-mass energies, y/s NN 200 GeV, available in 
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Au+Au collisions at the Relativistic Heavy Ion Collider 
(RHIC), have resulted in a significant production of high 
transverse momentum hadrons (px > 2 GeV/c) for the 
first time in heavy-ion physics. High pt particle produc- 
tion in hadronic collisions results from the fragmentation 
of quarks and gluons emerging from the initial high Q 2 
parton-parton scatterings [2]. Thus, hard processes in 
nucleus-nucleus (AA) collisions provide direct informa- 
tion on the early partonic phases of the reaction. In the 
absence of nuclear medium effects, hard scattering yields 
in AA reactions are expected to scale like an incoherent 
superposition of nucleon-nucleon (NN) collisions because 
of the small probability of hard scattering processes per 
NN collision. In AA reactions, the number of NN col- 
lisions (N co u) at impact parameter b is simply propor- 
tional to the geometric nuclear overlap function, TaaQ>), 
and can be calculated in an eikonal approach [3]. Af- 
ter scaling by the nuclear geometry, spectra of high pt 
particles measured in AA reactions can be compared to 
the baseline "vacuum" (p + p, p + p) and "cold medium" 
(p, d+A) data, as well as to perturbative [4-7] and "clas- 
sical" [8] QCD predictions. Any departure from the ex- 
pected 7V co ;/-scaled result provides information on the 
strongly interacting medium in central heavy-ion reac- 
tions. 

One of the most significant observations from the first 
RHIC run (Run-1) was the suppressed yield of moder- 
ately high pt neutral pions (pt ~ 1-5 - 4.0 GcV/c) in 
central Au+Au at y/s NN = 130 GeV with respect to 
the iV co ;/-scaled p + p and peripheral Au+Au data [9] . 
This result points to strong medium effects present in 
central Au+Au and has triggered extensive theoretical 
studies on its origin [4-8,10-13]. Most of these studies 
are based on the prediction [14,15] that a QGP would 
induce multiple gluon radiations from the scattered fast 
partons, effectively leading to a suppression of high pt 
hadronic fragmentation products ("jet quenching"). Al- 
ternative interpretations have been proposed based on 
initial-state gluon saturation [8] or final-state hadronic 
interactions [13]. 

This Letter presents ir° results obtained by the 
PHENIX experiment in Au+Au collisions at v /s NN = 
200 GeV and compares them to the p + p — > tt° X data 
measured in the same experiment at the same center- 
of-mass energy [16]. During the 2001-2002 RHIC run 
(Run-2), PHENIX obtained ir° data measured in the elec- 
tromagnetic calorimeter (EMCal). The present analysis 
uses 30 million minimum bias events, triggered by a co- 
incidence between the Zero Degree Calorimeters (ZDC) 
and the Beam-Beam Counters (BBC), with vertex posi- 
tion \z\ < 30 cm. In Run-2, the PHENIX EMCal was 
fully instrumented providing a total solid angle coverage 
at mid-rapidity of approximately A?? = 0.7 and A<f> = tt 
and the total collected n° statistics was a factor of ~ 100 
larger than in Run-1 [9]. The combination of larger ac- 
ceptance, high statistics, and the measurement of p + p 



data in the same detector permits a precise study of the 
high pt tt° production mechanisms in AA collisions at 
RHIC. 

Neutral pions are reconstructed via their 7r° — > 77 de- 
cay through an invariant mass analysis of 7 pairs detected 
in the EMCal [17] which consists of 6 lead-scintillator 
(PbSc) and 2 lead-glass Cerenkov (PbGl) sectors. The 
large radial distance of the calorimeters to the interaction 
region (> 5 m) and their fine granularity (A77 x A<f> 
0.01x0.01) keep the tower occupancy low, <10% even 
in the highest multiplicity Au+Au events. The energy 
calibration is obtained from beam tests and, in the case 
of PbSc, from cosmic rays data and the known mini- 
mum ionizing energy peak of charged hadrons travers- 
ing the calorimeter. It is then confirmed using the n° 
mass, as well as the agreement of the calorimeter en- 
ergy with the measured momentum of identified elec- 
trons. The systematic error on the absolute energy scale 
is less than 1.5%. Photon-like clusters are identified in 
the EMCal by applying time-of-flight and shower profile 
cuts [17]. The selected clusters are binned in pair invari- 
ant mass m 77 and pt- An additional energy asymmetry 
cut, \E-yt - E- y2 \/(E- yl + E l2 ) < 0.7(PbGl), 0.8(PbSc), 
is applied to the reconstructed pairs. The signal-to- 
background in peripheral (central) is approximately 20 
(5) and 0.5 (0.01) for the highest and lowest pt, respec- 
tively. The combinatorial background is estimated and 
subtracted by mixing clusters from different events with 
similar centrality and vertex, and normalizing the dis- 
tribution in a region outside the n° mass peak. The n° 
yield in each pt bin is determined by integrating the sub- 
tracted m 77 distribution in a ±3<t window determined by 
a pT-dependent parameterization of Gaussian fits to the 
7T peaks. 

The raw PbSc and PbGl 7r° spectra are normalized to 
one unit of rapidity and full azimuth (this acceptance 
correction rises quickly with pt up to a ~l/0.35 pure ge- 
ometric factor). The spectra are further corrected for (i) 
the detector response (energy resolution, dead areas), (ii) 
the reconstruction efficiency (analysis cuts), and (iii) the 
occupancy effects (cluster overlaps). These corrections 
are quantified by embedding simulated single 7r°'s from 
a full PHENIX GEANT [18] simulation into real events, 
and analyzing the merged events with the same analysis 
cuts used to obtain the real yields. Each correction is 
determined, for each centrality bin, as the ratio of the 
input (simulated) to the reconstructed pt distribution. 
The overall yield correction amounts to ~2.5 with a cen- 
trality dependence <25%. The losses are dominated by 
fiducial and asymmetry cuts. 

The main sources of systematic errors in the PbSc and 
PbGl measurements are due to the uncertainties in: (i) 
the yield extraction, (ii) the yield correction, and (iii) the 
energy scale. The relative contributions of these effects 
to the total error differ for the PbSc and PbGl (Table I). 
The weighted average of the two independent measure- 
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ments reduces the total error. The nominal energy reso- 
lution [17] is adjusted in the simulation to reproduce the 
true width of the tt° peak observed at each px, smear- 
ing the energies with a constant term of 7% for PbSc 
and —9% for PbGl. The shape, position, and width of 
the 7T° peak measured in all different centralities are then 
confirmed to be well reproduced by the embedded data. 
The final systematic errors on the spectra are at the level 
of -10% at 1 GcV/c and -17% at the highest p T (Ta- 
ble I). A correction for the true mean value of the pt 
bin is applied to the steeply falling spectra. No correc- 
tions have been applied to account for the contribution 
of feed-down 7r°'s (mainly coming from and r\ decays) 
which are <5% based on HIJING [19] simulations. 

The event centrality is determined by correlating the 
charge detected in the BBC with the energy measured in 
the ZDC detectors. A Glauber model Monte Carlo (MC) 
combined with a simulation of the BBC and ZDC re- 
sponses [20-22] gives an estimate of the associated num- 
ber of binary collisions (N co u) and participating nucleons 
(Np a rt) in each centrality bin (Table II). The fully cor- 
rected and combined PbSc and PbGl 7r° pt distributions 
are shown in Fig. 1 for minimum bias and for 9 centrality 
bins scaled by factors of 10 for clarity. 

We quantify the medium effects on high pt produc- 
tion in AA collisions with the nuclear modification factor 
given by the ratio of the measured AA invariant yields 
to the AW collision scaled p + p invariant yields: 



Raa{pi 



(l/N^)d 2 Nf A /d PT dy 
{N co u)/(J% el x cPog/dprdy' 



(1) 



where the (N co u) / 'a™ is just the average Glauber nu- 
clear overlap function, {TauAu), in the centrality bin un- 
der consideration (Table II). Raa(pt) measures the de- 
viation of AA data from an incoherent superposition of 
NN collisions. For px ^ 2 GeV/c, Raa is known to 
be below unity, since the bulk of particle production is 
due to soft processes which scale closer to the number of 
participant nucleons [20] than to N co u. 

Figure 2 shows Raa a s a function of px for ir° mea- 
sured in 0-10% central (closed circles) and 80-92% periph- 
eral (open circles) Au+Au. The PHENIX p + p — > 7r° 
data [16] is used as the reference in the denominator. 
The Raa values for central collisions are noticeably be- 
low unity, as found at 130 GeV [9], and in contrast to 
the enhanced high px 7r° production (Raa >1) observed 
at CERN-SPS energies [23] and interpreted in terms of 
initial-state px broadening effects ("Cronin effect" [24]). 
Within errors, peripheral Au+Au collisions behave like 
a superposition of p + p collisions with regard to high 
Pt 7t° production (Raa ~ !)■ In central collisions, the 
suppression is smallest at 2 GeV/c and increases to an ap- 
proximately constant suppression factor of 1/Raa ~ 4-5 
over the pt range of 4 - 10 GeV/c, —30% above the ex- 
pectation from N P art scaling (dotted line in Fig. 2). 



The magnitude and pt dependence of Raa (corre- 
sponding to parton fractional momenta x ~ 2px/\/s — 
0.02 — 0.1 at midrapidity) , is inconsistent with the ex- 
pectations of leading-twist "shadowing" effects on the 
nuclear parton distribution functions alone [25]. Dif- 
ferent jet quenching calculations [4-7,10-12], based on 
medium-induced radiative energy loss, can reproduce the 
magnitude of the ir° suppression assuming the forma- 
tion of a hot and dense partonic system. The predicted 
Pt dependence of the quenching, however, varies in the 
different models. All models that include the Landau- 
Pomeranchuck-Migdal (LPM) interference effect [15,26] 
predict Raa effectively oc y/pr [10]. Such a trend is not 
compatible with our data over the entire px range. Anal- 
yses which combine LPM jet quenching together with 
shadowing and initial-state px broadening generally re- 
produce the whole pt dependence of the ir° suppres- 
sion [4], as do recent approaches that take into account 
detailed balance between parton emission and absorp- 
tion [7]. However, based solely on the data presented 
here, we are not able to distinguish between partonic or 
hadronic [13] energy loss scenarios. 

The centrality dependence of the high pt tt° suppres- 
sion is shown in Fig. 3 as a function of N part . The 
suppression is characterized as the ratio of Au+Au over 
p + p yields integrated above 4 GeV/c and normalized 
using two different scalings. Raa (circles) denotes the 
N co u scaling as in Eq. (1), whereas R V aa (crosses) indi- 
cates Np ar t scaling expected in scenarios dominated ei- 
ther by gluon saturation [8] or by surface emission of the 
quenched jets [10]. Fig. 3 indicates that the transition 
from the N co u scaling behaviour (Raa — 1) apparent in 
the most peripheral region, to the strong suppression seen 
in central reactions (Raa ~ 0.25) is smooth. In addition, 
although there is no exact participant scaling (R V aa > 
1 for all centralities), the ir° production per participant 
pair above 4 GeV/ c is approximately constant over a wide 
range of intermediate centralities, in qualitative agree- 
ment with a parton saturation model prediction [8] . 

In summary, transverse momentum spectra of neutral 
pions have been measured at mid-rapidity up to pt ~ 
10 GeV/c for 9 centrality bins of Au+Au collisions at 
y/s NN = 200 GeV. The spectral shape and invariant yield 
for peripheral reactions are consistent with those of p + p 
reactions scaled by the average number of inelastic NN 
collisions. Central yields, on the other hand, are signif- 
icantly lower than peripheral Au+Au and p + p scaled 
yields, as found at t/s^ = 130 GeV. The observed sup- 
pression increases slowly with pt to as much as a factor 
4-5 in the 10% most central collisions, remaining con- 
stant within errors above —4 GeV/c. The suppression 
is already apparent in semi-peripheral reactions and in- 
creases smoothly with centrality. The magnitude of the 
deficit can be reproduced by parton energy loss calcu- 
lations in an opaque medium, but its pt and centrality 
dependence puts strong constraints on the details of en- 
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ergy loss and the properties of the medium. The role of 
initial-state effects, including shadowing, px broadening, 
and gluon saturation will be studied with data from the 
recent RHIC run using d+Au, where final-state medium 
effects such as jet quenching are absent. 
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TABLE I. Summary of the dominant sources of systematic 
errors on the PbSc and PbGl tv° yields, and total errors on 
the combined measurements. The error ranges are quoted 
for the lowest to highest pr values. For the combined it 
spectra and Raa, we quote the approximate i) statistical, ii) 
point-to-point systematic, and iii) absolute normalization (for 
most central and peripheral reactions) errors. 



Source 




syst. error PbSc 


syst. error PbGl 


Yield extraction 




10% 


6 - 7% 


Yield correction 




8% 


8% 


Energy scale 




3 - 11% 


7 - 13% 


Total error (%) 


stat. 


syst. 


normalization 








central peripheral 


Spectra 


2 - 40 


10 - 17 


5 5 


Raa 


2 - 45 


11 - 22 


14 30 



TABLE II. Centrality bin, average nuclear overlap func- 
tion, number of NN collisions, and number of participant 
nucleons obtained from a Glauber MC [21,22] and the BBC 
and ZDC responses for Au+Au at ^/snn = 200 GeV. The 
centrality bin is expressed as percentiles of oauAu = 6.9 b. 
The last line refers to minimum bias collisions. 



Centrality 


(TauAu) (mb" 1 ) 


(N coll ) 


(Apart) 


0-10% 


22.75 ± 1.56 


955.4 ± 93.6 


325.2 ± 3.3 


10-20% 


14.35 ± 1.00 


602.6 ± 59.3 


234.6 ± 4.7 


20-30% 


8.90 ± 0.72 


373.8 ± 39.6 


166.6 ± 5.4 


30-40% 


5.23 ± 0.44 


219.8 ± 22.6 


114.2 ± 4.4 


40-50% 


2.86 ± 0.28 


120.3 ± 13.7 


74.4 ± 3.8 


50-60% 


1.45 ± 0.23 


61.0 ± 9.9 


45.5 ± 3.3 


60-70% 


0.68 ± 0.18 


28.5 ± 7.6 


25.7 ± 3.8 


60-80% 


0.49 ± 0.14 


20.4 ± 5.9 


19.5 ± 3.3 


60-92% 


0.35 ± 0.10 


14.5 ± 4.0 


14.5 ± 2.5 


70-80% 


0.30 ± 0.10 


12.4 ± 4.2 


13.4 ± 3.0 


70-92% 


0.20 ± 0.06 


8.3 ± 2.4 


9.5 ± 1.9 


80-92% 


0.12 ± 0.03 


4.9 ± 1.2 


6.3 ± 1.2 


min. bias 


6.14 ± 0.45 


257.8 ± 25.4 


109.1 ± 4.1 
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FIG. 1. Invariant 7r° yields at mid-rapidity as a function of 
pr for minimum bias and 9 centralities in Au+Au at y/s NN 
= 200 GeV [0-10% (80-92%) is most central (peripheral)]. 
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FIG. 2. Nuclear modification factor Raa(pt) for ty° in cen- 
tral (closed circles) and peripheral (open circles) Au+Au at 
^s NN = 200 GeV. The error bars include all point-to-point 
experimental (p + p, Au+Au) errors. The shaded bands rep- 
resent the fractional uncertainties in (TauAu) and in the tv° 
yields normalization added in quadrature, which can move all 
the points up or down together (in the central case the shaded 
band shown is the fractional error for the first point). 
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FIG. 3. Ratio of Au+Au over p + p ir° yields integrated 
above 4 GeV/c and normalized using N co u (circles) and N par t 
(crosses), as a function of centrality given by N par t- The er- 
rors bands and bars are the same as for Fig. 2. 
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